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Abstract 

Calculations  of  the  effect  of  bulk  traps  on  the  time  dependence  of  the 
channel  current  in  InP  accumulation-type  MISFET's  are  reported.  It  is  shown 
that  a  uniform  density  of  deep  acceptors  in  the  bulk  material  given  rise  to  a 
decay  in  the  channel  current  which  decays  approximately  linearly  as  a  function 
of  the  logarithm  of  time.  If  capture  of  the  electrons  is  assumed  to  occur  by 
a  nonradiative  multi  phonon  emission  process,  then  the  temperature  dependence 
that  is  obtained  agrees  with  experimental  results. 
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Introduction 


Indium  phosphide  (InP)  is  an  attractive  material  for  the  fabrication  of 
metal -insulator-semi conductor  (MIS)  field-effect  transistor  (FET)  because  of 
its  large  saturated  drift  velocity  for  electrons  and  the  favorable  conditions 
present  at  its  oxide-semiconductor  interface.  While  efforts  to  construct 
MISFETs  on  gallium  arsenide  were  unsuccessful  due  to  the  large  interface  state 
density  on  this  material^  pr0gress  has  been  much  more  rapid  in  the  case 
of  In?.  Interface  state  densities  between  deposited  layers  of  silicon  dioxide 
(SiC^)  and  InP  have  been  calculated  to  be  as  low  as  2  x  10*1  cm“^eV*^  based  on 
capacitance-voltage  measurements  on  MIS  diodes^>4»^.  A  type  of  MISFET  which 
has  been  labeled  the  accumulation-type  device  has  been  realized  on  InP^J. 
Construction  of  this  device  begins  with  the  formation  of  two  heavily  doped 
(n+)  ohmic  contacts  for  the  source  and  drain  regions  on  a  wafer  of  semi- 
insulating  (SI)  InP.  The  wafer  is  then  coated  with  a  thin  insulating  layer* 
overlapping  gate  electrodes  are  evaporated  on  top  of  the  insulator  and  windows 
are  etched  through  the  insulator  for  connection  to  the  source-drain 
electrodes.  A  cross-sectional  view  of  the  device  is  shown  in  Fig.  1.  This 
device  was  originally  proposed  and  fabricated  oecause  Che  reduced  capacitance 
between  the  conducting  channel  and  the  substrate  was  expected  to  lead  to  a 
better  high  frequency  response  than  similar  devices  constructed  on  p-type 
substrate  material.  Indeed,  this  was  found  to  be  the  case. 

A  remaining  problem  to  be  solved  before  the  device  is  perfected  is  a 
long-term  drift  in  the  drain  current  that  is  observed  when  a  constant  voltage 
is  applied  between  the  gate  and  source  electrodes^ Several  models  have 
been  proposed  which  attempt  to  explain  the  origin  of  the  drift  by  tunneling 
currents  between  states  in  the  insulator  and  the  conduction  band  of  the  semi - 
conductor^  j  feel  that  these  models  have  serious  shortcomings  and 
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Abstract 


The’ enclosed  reports  represent  work  performed  at  UCSD  on  Contract 

- «-r.  i 

N00014-82-K-2032  entitled  ‘Surface  and  Interfacial  Properties  of  InP*  and 
provides  a  full  account  of  the  results  obtained  during  the  contract  period: 

Hay  1,  1984  through  April  31,  1985.  The  paper  ‘Space  charge-limited  currents 
and  trapping  in  semi-insulating  InP"'  has  now  been  published  in  Electron. 

Device  Letters,  volume  EDL-6,  page  356  {1985).  The  manuscript  “Effect  of  bulk 
traps  on  the  In?  accumulation  type  MISFE!*1  will  be  presented  as  an  invited 
talk  at  the  fall  meeting  of  the  Electrochemical  Society  in  Las  Vegas,  October 
14-17  and  will  be  published  in  the  Journal  of  the  Electrochemical  Society. 


that  although  they  give  substantially  correct  mathematical  expressions  for  the 
current  drift  they  have  errors  in  their  formulation.  This  is  further 
discussed  in  the  paper. 

A  first  order  model  for  the  I-V  behavior  of  the  accumulation-type 
transistor  based  on  the  gradual  channel  approximation  has  been  formulated  by 
Wieder^^.  His  model  assumed  that  the  substrate  (the  SI  InP)  could  be 
treated  as  an  ideal  insulator. 

In  this  paper  I  present  calculations  for  the  time-dependent  behavior  of 
the  accumulation-type  transistor  which  take  into  account  the  fact  that  the 
conducting  channel  of  the  device  is  on  the  surface  of  a  semiconductor  which 
contains  both  shallow  donors  (due  to  native  defects)  and  deep  acceptors 
(introduced  by  iron  doping).  The  time  dependence  of  the  device  is  calculated 
by  assuming  that  the  capture  and  emission  of  electrons  by  the  bulk  impurity 
states  is  governed  by  a  Shockley-Read  type  of  recombi nation  process. 

Poisson's  equation  has  been  solved  numerically  to  obtain  the  time  development. 


Background 


The  key  features  which  must  be  explained  by  any  proposed  model  for  drift 
of  the  channel  current  in  the  accumulation-type  transistor  are  illustrated  in 
data  published  by  Lile,  et  al  .^3  which  are  reproduced  in  Fig.  2.  The  channel 
current  decays  linearly  as  a  function  of  the  logarithm  of  time  and  the  capture 
process  "freezes  out"  at  low  temperature.  Clearly  a  tunneling  process  as 
decribed  by  Heiman  and  Warfield^1^  cannot  explain  such  data  since  there  would 
be  no  temperature  dependence  if  electrons  were  to  tunnel  from  conduction  band 
states  to  oxide  trap  levels  at  the  same  energy. 

One  of  the  proposed  models  for  explaining  the  current  drift  in  the 
accumulation-type  transistor  as  presented  by  Goodnick,  et  al  attractive 
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because  it  is  based  on  physical  evidence  that  the  interface  of  InP  coated  with 

deposited  SiC^  consists  of  multiple  layers  as  shown  in  Fig.  3a.  They  proposed 

that  the  layer  immediately  adjacent  to  the  InP  consists  of  a  relatively  large 

bandgap  material  (InPO^)  which  is  covered  with  a  lower  bandgap  material 

(^203)  which  in  turn  is  covered  with  Si Og -  When  a  positive  voltage  is 

applied  to  the  gate  of  the  device  to  form  the  conducting  channel  in  the  SI 

InP,  they  proposed  that  a  thermally-assisted  tunneling  current  flows  from  the 

semiconductor  to  states  in  the  In203.  Their  solution  for  the  channel  current 

(Id)  as  a  function  of  time  gives  a  linear  dependence  of  Id  on  log  t.  The 

temperature  dependence  is  obtained  in  their  model  from  the  experimental ly 

measured  conduction  band  offset  between  the  InP  and  the  Wager  et 

al  .E1^]  conciuc|e(j  from  photoemission  (XPS,  UPS)  and  electron-loss  spectroscopy 

(ELS)  measurements  of  deposited  Si O2  layers  on  InP  that  the  conduction  band 

minimum  in  the  native  oxide  is  0.1  to  0.2  eV  higher  than  in  InP.  In  the 

Goodnick,  et  al .  thermally-assisted  tunneling  model,  the  electrons  can  tunnel 

from  the  InP  to  the  ^03  only  if  their  kinetic  energy  is  greater  than  the 

energy  barrier  step  from  the  InP  to  the  ^03.  As  the  temperature  is  reduced 

fewer  electrons  are  able  to  surmount  this  barrier  and  the  drift  effect 

decreases.  However,  this  model  fails  to  take  into  account  the  fact  that  with 

a  positive  voltage  on  the  gate  electrode  there  will  be  a  voltage  drop  across 

the  InPO/j  layer  and  a  lowering  of  the  conduction  band  of  the  10203  as  shown  in 

Fig.  3b.  Once  the  conduction  band  of  the  ^03  is  below  the  conduction  band 

of  the  InP  all  the  temperature  dependence  of  the  Id  vs.  t  curves  would 

disappear.  For  a  40  A  layer  of  InPQ^  as  assumed  in  their  calculation  this  requires 

a  gate  voltage  of  -  .  2(d^.Q  +  ^InPO  ^rfInP0  =  ^  volts*  1°  This  very  rough 

2  4  4 

estimate,  d^  represents  the  thickness  of  the  InPO^  and  Si02  layers,  the  dielectric 
constants  of  the  two  layers  have  been  assumed  equal,  and  the  amount  of  conduction 
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band  offset  has  been  assumed  to  be,  as  an  upper  limit,  0.2  v.  The  temperature 
dependence  of  the  drift  current  should  then  disappear  for  gate  voltages  larger  than 
about  5  v.  This  is  contrary  to  the  experimental  evidence.  At  low  temperature  the 
current  drift  disappears  for  all  values  of  gate  voltage.  One  could  attempt  to 
salvage  this  model  by  reducing  the  thickness  of  the  InPO^  layer  so  as  to 
minimize  the  voltage  drop  across  it.  However,  this  would  at  the  same  time 
eliminate  the  tunneling  barrier  and  eliminate  the  log  t  dependence  of  1^. 
Alternatively,  one  could  assume  a  lower  resistivity  of  the  InP04  layer  so  that 
the  voltage  drop  across  it  would  be  less,  but  this  would  again  eliminate  the 
tunneling  behavior  needed  to  get  the  proper  time  dependence. 

Okamura  and  Kobayashi^1^  have  proposed  that  the  current  drift  in  the 
accumulation  type  device  can  be  explained  by  thermally-assisted  tunneling  from 
the  conduction  band  of  the  InP  into  a  pair  of  electron  trap  levels  in  the 
insulator.  The  energy  band  diagram  relevant  to  their  scheme  is  shown  in 
Fig.  4.  They  base  their  calculation  on  a  model  originally  proposed  by 
Koelmans  and  DeGraafft15^.  This  model  begins  by  assuming  that  the  charge  flow 
into  the  oxide  trap  can  be  represented  by  a  Shockley-Read^16^  type  of 
recombination  rate  equation,  viz: 

dntr 

ST  *  s(x)  v  ["s  (Ntr  ’  ntr>  '  "l  "tr1  (1> 


where  n$  is  the  electron  concentration  at  the  interface,  0|.r  is  the  density  of 
filled  traps,  V  is  the  average  thermal  velocity  of  the  electrons  and  Ntr  is 
the  total  trap  density.  The  quantity  nj ,  is  given  by  the  equation 


Nc  e-(Ecs  -  E(x))/kT 


(2) 
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where  Nc  is  the  effective  density  of  states  for  the  conduction  band,  k  is 
Boltzmann's  constant,  T  is  the  absolute  temperature,  and  Ecs  and  E(x)  are 
energies  defined  in  Fig.  4.  This  particular  value  of  nj  must  be  used  if  the 
principle  of  detailed  balance  is  to  be  maintained.  That  is,  at  thermal 


equi li bri urn 


(3) 


which  then  implies  that 


Si  nee 


and 


ns  (Ntr  "  "tr> 


ns  =  Nc  e 


(ECS)AT 


Ntr  ~  ntr  eE(x)/kT 
ntr 


(4) 

(5) 

(6) 


analagous  to  the  equation  for  the  ratio  of  ionized  donors  to  unionized  donors 
for  an  impurity  state  in  a  semiconductor^-17-!,  one  is  forced  to  choose  Eq.  (2) 
to  represent  nj  in  order  to  satisfy  the  detailed  balance  requirement.  The 
tunneling  current  is  modeled  by  assuming  that  the  charge  flow  into  the  oxide 
traps  can  be  accounted  for  by  a  capture  cross  section  which  is  an  exponential 
function  of  distance  into  the  oxide,  i.e. 


S(x)  =  S0  e-ax  (7) 

where  1/a  represents  a  tunneling  length.  In  this  picture  the  probability  that 
an  electron  at  the  surface  will  tunnel  into  an  oxide  trap  is  represented  by 
the  product  of  the  electron  density  at  the  surface,  the  thermal  velocity,  the 
number  of  empty  traps  at  distance  x  into  the  oxide,  and  the  distance-dependent 
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capture  cross  section.  Koelmans  and  OeGraaff  solve  this  equation  to  obtain  an 

equation  for  the  total  trapped  charge  density,  (N  (t),  as  a  function  of  time 

e 

(t)  which  is  of  the  form 


N  (t)  =  —  In  (t/t  ) 
e  '  a  o 


where  tQ  is  a  constant.  The  original  model  as  formulated  by  Koelmans  and 

DeGraaff  gives  no  temperature  dependence  to  the  tunneling  current. 

The  modification  to  this  model  as  proposed  by  Okamura  and  Kobayashi 

includes  a  trap  level  in  the  oxide  which  is  placed  above  the  conduction  band 

(2) 

of  the  InP  by  an  energy  (E§v  '  -  E  ).  They  reasoned  that  only  the  more 
energetic  conduction  band  electrons  would  be  able  to  be  captured  by  these 
traps  and  that  the  capture  rate  would  be  reduced  by  a  factor 
e~^s^  “  Ecs>/kT.  If  the  Shockley-Read  equation  is  solved  using  their 
capture  cross  section,  then  the  charge  trapped  by  this  level  is  given  by 

N  <2>  -E  (2)AT 

Ne-V  S  l»(t/t#).  (9) 
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However,  their  derivation  neglects  the  constraint  that  detailed  balance 
imposes  on  the  formulation  of  the  Shockley-Read  rate  equation.  If  the  rate 
equation  is  required  to  satisfy  the  requi rement  that 


(2) 

tr _  (10) 

dt  "  u 


(2) 

at  thermal  equilibrium  then  the  term 
must  be  changed  to 


in  the  emission  rate  expression 


7 


"l(2)  -  Nc  e 


-<V  '  -  Ecs)AT  -(Ecs  -  E(x))AT 


1  N„  e 
c 


-(Es<2)  -  e(x))/kT 


(U) 


This  is  mathematically  equivalent  to  assuming  a  smaller  capture  cross  section 
for  the  trap.  When  the  rate  equation  is  solved  with  this  reformulation  there 
is  again  no  temperature  dependence  of  the  trapped  charge  density.  Thi sis 
consistent  with  the  paper  of  Koelmans  and  DeGraaf since  they  formulated 
their  model  assuming  that  the  trap  level  was  placed  at  an  arbitrary  position 
in  the  forbidden  energy  gap  and  obtain  the  generalized  result,  that  the 
trapping  of  electrons  in  the  oxide  is  not  temperature  dependent.  The  question 
of  whether  or  not  one  is  justified  in  neglecting  the  detailed  balance 
requirement  will  be  addressed  later  in  the  paper. 

Yamaguchi ,  et  al  have  calculated  the  tunneling  current  into  the 

oxide  using  a  model  similar  to  the  original  Koelmans  and  DeGraff  model. 
Although  they  employ  a  different  calculation  scheme,  their  model  like  the 
Koelmans  and  DeGraff  model  contains  no  element  for  predicting  the 
experimentally  observed  temperature  dependence  of  the  drift  current. 

Given  the  problems  that  have  been  encountered  in  devising  a  model  which 
simultaneously  gives  the  correct  time  and  temperature  of  the  drift  current  in 
the  accumulation-type  transistor  when  the  charge  trapping  mechanism  is  placed 
in  the  insulator,  it  is  interesting  to  consider  an  alternate  approach  which 
places  the  charge  trapping  mechanism  in  the  semiconductor.  How  do  we  know 
that  the  charge  is  being  trapped  in  the  insulator?  The  answer  is  that  we 
don '  t .  The  extant  published  data  do  not  rule  out  the  possibility  that  the 
charge  trapping  may  occur  in  the  semiconductor.  The  next  point  to  consider  is 
that  the  accumulation-type  device  has  as  its  basis  an  electric  field  induced 
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layer  of  electrons  on  the  surface  of  an  InP  SI  crystal.  The  InP  substrate 
material  has  a  background  concentration  of  unintentionally  introduced 
hydrogenic  donors  which  are  located  7  to  10  meV  below  the  conduction  band 
minimum  (£c)  [18]  and  have  a  density  which  is  of  the  order  of  2x10  cm-.  . 

In  the  near  surface  region,  phosphorus  vacancies  may  also  exist.  These  are 
thought  to  give  rise  to  an  energy  level  0.1  eV  below  Ec  [19].  In  these 
calculations  was  placed  0.1  eV  below  Ec  although  the  results  were  not 
particularly  sensitive  to  the  value  used  for  this  parameter.  The  origin  of 
these  native  donors  is  at  the  present  time  unknown.  During  the  growth  of  the 
crystal,  iron  is  added  to  the  melt  to  produce  deep  compensating  acceptors 
which  are  located  approximately  0.68  eV  above  the  valence  band  maximum 
(Ey).  As  long  as  the  deep  acceptor  density  is  greater  than  the  background 
concent  rati  on  of  native  donors,  the  material  will  be  semi -i nsulati ng.  As  the 
present  time,  commercially  available  SI  InP  is  not  sold  with  clearly  specified 
upper  limits  to  the  amount  of  iron  which  may  be  present.  Consider  now  Fig.  5 
which  illustrates  the  energy  band  diagram  of  an  SI  InP  substrate  biased  such 
that  the  surface  is  accumulated.  If  the  material  is  fully  compensated  and  the 
banJs  are  flat,  then  at  thermal  equilibrium  nearly  ail  the  donors  vN^  =  Nd) 
will  be  ionized  but  only  a  fraction  (N_^  *  N^)  of  the  acceptors 
will  be  ionized.  The  ionized  donor  density  will  be  given  by 


the  ionized  acceptor  density  will  be  given  by 
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the  electron  density  will  be  given  by 


-  EF)/kT 


and  the  hole  density  will  be  given  by 


”(Ep  -  E  J/kT 

P  -  Ny  e 


(14) 


(15) 


In  these  equations  g,  represents  the  degeneracy  of  states  for  donors  and 

0 

acceptors  respectively,  Ep  represents  the  position  of  the  Fermi  level  with 
respect  to  E  in  the  neutral  bulk  material,  E.  represents  the  energy  of  the 
donor  and  acceptor  levels  respectively  as  measured  with  respect  to  Ev,  Nc 
represents  the  effective  density  of  states  at  E  =  Ec  for  electrons  and  Ny 
represents  the  effective  density  of  states  at  E  =  Ey  for  the  holes.  The 
position  of  the  Fermi  level  can  be  readily  calculated  by  numerically  solving 
the  equation 


P  =  -n  +  P  +  Ndi  -  Ngi  =  0.  (16) 

for  Ep  and  using  equations  12,  13,  14  and  15  to  solve  for  the  desired 
quantities.  If  Na  is  greater  than  several  times  Nd  and  Ed  is  greater  than  Ey 
by  several  kT,  then  to  a  good  approximation 


E 


F 


(17) 


When  a  positive  voltage  is  applied  to  the  gate  electrode  of  an  MOS  device 


on  SI  InP,  the  energy  bands  will  bend  downward  as  shown  in  Fig.  5.  Notice 


that  acceptor  states  near  the  surface  that  were  formerly  above  the  Fermi  level 
are  now  pushed  downward  by  electric  field.  As  this  happens  electrons  will 
fill  the  empty  acceptor  states  causing  them  to  be  completely  ionized  in  the 
near  surface  region.  If  the  gate  voltage  is  held  constant,  then  the  electric 
field  which  now  terminates  on  the  additionally  ionized  acceptors  will  clearly 
cause  a  reduction  in  the  mobile  channel  charge.  This  process  would  lead  to  a 
reduction  of  the  channel  current  of  an  FET  over  time  and  could  easily  be 
mistaken  for  interfacial  charging.  The  remainder  of  this  paper  presents 
numerical  calculations  of  the  time  dependence  of  this  drift  effect  under 
various  conditions. 

Theory 

The  first  step  in  the  formulation  of  this  model  is  the  assumption  that  n+ 
ohmic  contacts  are  present  on  the  surface  of  the  semiconductor  and  that  the 
gate  electrode  partially  overlaps  these  contacts.  This  is  the  typical  state 
of  affairs  for  the  accumulation-type  MISFET.  This  assumption  is  essential 
from  a  physical  standpoint  because  the  n+  contact  provides  for  easy  transfer 
of  electrons  in  and  out  of  conduction  band  states  and  is  essential  for  the 
formulation  of  the  mdel  because  if  permits  th»  electron  density  to  be 
described  by  an  equilibrium  Fermi  level.  The  generation  recombi nation 
processes  were  assumed  to  be  characterized  by  the  Shockley-Read  rate 
equations^6^. 


dN.,  (x)  mb/kT  Nr  "  Er)/kT 

-& - Cd[-»ljW.b^b{e  -*„(x)i3 


'di 


(13) 


and 


d_!i^Lc  [|N  -N  e**1 

dt  a  a  ai  b 


N 


c  (Ea  ‘  Ec)/kT 
-e  a  C  N  .(x)] 

2  a1 


(19) 
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where  and  Ca  represent  capture  coefficients  given  by  vcj^  and  voa 
respecti  vely ,  n^  represents  the  electron  density  in  the  bulk  material,  <{>  is 
the  static  potential  and  a ^  and  oQ  represents  the  capture  cross  sections  for 
the  donors  and  acceptors  respectively.  Both  the  donor  and  acceptor  impurity 
states  were  assumed  to  be  doubly  degenerate.  The  boundary  conditions  were  as 
follows:  it  is  assumed  that  the  fixed  interfacial  charge  is  equal  to  zero  and 
that  the  work  function  difference  between  the  metal  and  the  semiconductor  is 
equal  to  zero.  A  gate  voltage  Vg  =  0  is  applied  for  t<0  and  a  fixed  positive 
gate  voltage  is  applied  for  t  >  -0.  With  these  boundary  conditions 

NHi (x)  j  =  N  and  N.(x)|  =  N..  (20) 

Qi  lt=g  a  ai  ltsig  a 


The  value  of  the  surface  potential  at  t  =  0+  was  obtained  by  solving  the 


equati on 


Q  +  0  e  E  0 

V  =  *s  +  - —  =  *s  +  -I*1-  +  (T1 

ox  ox  ox 


where  Qs  and  Es  are  respectively  the  total  charge  on  the  semiconductor  and  the 

surface  electric  field,  Qss  is  the  charge  in  surface  states  and  CQX  is  the 

oxide  capacitance  per  unit  area.  For  the  calculated  results  presented  herein 

it  his  been  assumec.  that  Q$s  =  0.  The  effect  of  a  non-zero  threshold  voltage 

(Vt^)  can  be  obtained  in  this  model  by  assuming  a  non-zero  fixed  value  of 

Q^_.  The  boundary  condition  Vn  =  constant  would  then  be  replaced  with  the 

q  y 

boundary  condition  V  -  =  constant.  At  t  =  0+  the  ionized  donor  and 

^  'ox 

acceptor  densities  were  assumed  not  to  have  changed  from  their  values  for  t  < 
0.  The  equation  for  the  surface  electric  field  as  shown  below^-2^ 


12 


(-  q*  +  Ec  -  Ed  -  Ep)/kT 


+  Ndln  [?  +_?.(  y/kT. 

2  +  e  c  0  T 


] 


(Ef  -  Ea  -  Ey  4  q*)/kT  j 

+  Naln  I  —  ---E—E-v)-/kT  ~  ]  +  pb^e  '  q^AT- 


(22) 


2  +  e 


reduces  to 


E  =  -  ^1=  ( _ 

dx  1  J 


2kT  n 


lh 


aq*/2kT 


(23) 


for  q*  »  kT ;  pb>  nb  «  Nfl,  Nrf ,  and  Ndi  -  Nai . 

The  charge  density  distribution  at  t  =  0+  can  be  obtained  explicitly  by 
simultaneously  solving  Eqs.  21  and  23.  Similarly  the  charge  density 
distribution  at  t  =  “  can  be  obtained  by  simultaneously  solving  Eqs.  21  and 
22.  A  graphical  solution  of  these  equations  is  indicated  in  Figs.  6  and  7. 

The  Es  vs  <J»s  curves  at  t  =  0+  and  t  =  “  were  calculated  for  specific  donor  and 
acceptor  impurity  densities.  The  intersection  of  these  curves  with  a 
particular  curve  of  constant  Vg  as  formulated  by  Eq.  21  represents  a 
particular  solution.  The  area  between  the  Es  vs  i|>s  curves  at  t  =  0+  and  t  =  ® 
represents  the  field  of  possible  solutions  for  the  time  dependent  problem. 

The  decrease  in  the  conducting  channel  carrier  density  can  be  obtained  by 
noting  the  solution  for  ^  as  indicated  in  Figs.  6  and  7.  At  t  *  0+  the  total 
charge  due  to  conduction  band  electrons  is  given  to  a  good  approximation  by 

1,  q*/2kT 

Q_  *  e  E  *  (2kTe  nj*  e  S  .  (24) 

n  s  s  $  b 


At  t  =  “  the  mobile  electron  charge  can  be  obtained  by  numerically  integrating 
the  equation 

CO 

Qn  =  q /  n(x)dx 


(25) 


or  its  equi valent 


(26) 


where  ci4»  /dx  is  given  by  Eq.  22  and  is  obtained  from  from  Figs.  6  or  7. 
Although  Fig.  6  seemingly  indicates  that  the  solutions  differ  very  little  for 
large  Vg,  this  is  not  quite  the  case.  The  solutions  for  i>s  do  differ  by  a 
small  amount  but  since  the  mobile  electron  density  depends  exponentially 
on  '{'s,  the  difference  in  electron  density  can  be  large.  The  reduction  in 
channel  current  for  a  range  of  acceptor  concentrations  is  indicated  in  Fig.  8. 
It  is  instructive  to  consider  calculated  plots  of  the  surface  potential  and 
electric  field  in  Fig.  9  as  a  function  of  distance  into  the  crystal.  An 
unobvious  point  is  that  the  electric  field  penetrates  an  extremely  large 
distance  (w  -  2mm)  into  the  material  when  the  gate  voltage  is  initially 
applied.  This  is  due  to  the  fact  that  initially  only  a  small  fraction  of  the 
acceptors  are  ionized  and  there  are  few  charged  centers  for  the  electric  field 
lines  to  terminate  on.  If  a  small  geometry  FET  were  being  considered,  the 
entire  surrounding  area  of  the  FET  would  be  affected  by  the  field.  Only  after 
the  deep  acceptor  levels  near  the  surface  have  filled  with  electrons  does  the 
accumulation  layer  width  shrink  to  device  dimensions.  The  calculations 
presented  here  being  one-dimensional  in  nature  do  not  take  into  account  such 
two  dimensional  geometry  effects  since  the  InP  crystal  has  been  assumed  to  be 
semi  -i  nfi  ni  te. 

Solution  of  the  time  dependence  of  this  problem  begins  by  evaluating  Eq. 
22  at  ^  t  ,  inserting  this  into  Eq.  21  and  solving  the  result  iteratively 
using  Newton's  method  to  obtain  >1^.  The  variation  of  with  position  at  t  = 

0+  was  obtained  by  integrating  Eq.  23 
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/  „  ,  S  -q*kT 

/  dx  =  "  (jjFTrT'  '  e 
0  "  b  <|> 


d<l> , 


(27) 


which  yields 


2e  kT  !/>  .  ,  x  /ol  T  -q<f>  /2kT 

=  (-*  )^(e-t’*(x>/2kT  -e  5  ) 


"h<l 


(28) 


which  can  be  solved  for  ^(x )  to  give 


♦  (*,  =  (29) 

The  carrier  density  at  t  =  0+  is  then  obtained  by  inserting  Eq»  29  into 


n(x)  =  nb  e 


q’Hx)/kT 


(30) 


The  above  equations  were  used  to  solve  for  the  values  of  4>(x)  and  n(x)  on  an 

array  of  points  x^  with  a  variable  spacing  selected  with  the  algorithm 
-1°  _io 

xo  =  0,  xi  =  10  m,  Axx  =  10  m,  x^+1  =  x^  +  Ax;  Ax^+J  =  Ax;  if 
n (xi +1 )  *  l-ln(x^)  and  Ax^+1  =  l.lAx;  if  n(x..+1)  <  1.1  n(x).  The  variable 
spacing  is  required  in  order  to  have  closely  spaced  points  near  the  surface 
where  the  carrier  density  drops  rapidly  and  widely  spaced  points  in  the 
interior  where  the  carrier  density  drops  slowly. 

Given  initial  values  for  n(x),  Nd^(x)  and  Na^(x),  the  time  development  of 
the  trapped  charge  was  approximated  by  using  Eqs.  18  and  19  to  obtain 

dN  <x,t) 

Ndi(x,t  +  fit)  -  Ndi(x,t)  +  - fit  (31) 

and 

dN  .  (x,t) 

Nai(x,t  +  fit)  -  Nai(x,t)  +  - fit.  (32) 
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With  these  new  values  of  ionized  impurity  density  Poisson’s  equation 


*4--  P/s.  (33) 

dx  5 

can  be  solved  numerically  using  a  fourth-order  Runge-Kutta  met  hod  given 
that  the  charge  density  is 


P  =  q(-  n  +  Nd.(x,t  +  St)  -  Nai(x,t  +  St)) 

where 

n  =  eq*(x)/kT. 

The  second-order  differential  equation  can  be  reduced  to  the  pair  of  coupled 
fi  rst-order  equations  given  by 


(34) 

(35) 


and 


35 r-'V  (36) 


In  the  solution  of  Eq.  32  for  a  new  value  of  time,  say  t  =  t  +  St,  the 
first  trial  value  of  ♦  (t  +  St)  was  set  to  ts(t).  This  value  of  was  used 
in  Eq.  20  to  solve  for  the  trial  value  of  Es(t  +  St).  These  starting  values 
of  <J»s  and  E$  were  used  in  the  Runge-Kutta  procedure  to  solve  for  >k(x,t  +  St)  and 
E(x,t  +  St)  from  the  surface  inwards.  If  4* (x , t  +  St)  and  E(x,t  +  St)  reached 
zero  simultaneously  then  <|>  and  E-  represent  a  solution.  If  <f>(x,t  +  St) 
reached  zero  before  E(x,t  +  St),  then  ^  was  too  small  and  if  E(x,t  +  St) 
reached  zero  before  <J»(x,t  +  St)  then  <l»s  was  too  large.  A  new  value  of  ^  was 
chosen  that  was  known  to  lie  between  values  of  that  had  been  calculated  to 
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be  either  too  large  or  too  small  and  the  process  was  repeated  until  a  solution 
was  obtained. 

Calculated  results  for  the  charge  in  conduction  band  states  (Ns)  as  a 
function  of  time  are  shown  in  Fig.  10.  The  charge  is  nearly  constant  for  very 
short  times  and  then  exhibits  a  decay  which  is  nearly  linear  in  log  t.  The 
time  dependence  of  the  calculated  results  is  very  similar  to  the  data 
presented  by  Yamaguchi ,  et_alj8]  an(j  ^  Li  1  e *  et  al  .^1  However,  the 
temperature  dependence  is  not  as  strong  as  observed  in  the  measurements  by 
Taylor,  et  al .  The  temperature  dependence  in  the  model  just  described  is  due 
to  the  reduction  in  the  Debye  length 

2kTe 

i  _  r  S 


as  the  temperature  is  lowered.  At  low  temperature  the  electrons  crowd  closer 
together  at  the  surface  and  the  electron  density  in  the  far  surface  region  is 
reduced.  This  has  the  effect  of  reducing  the  magnitude  of  the  terms  due  to 
electron  capture  in  Eqs.  17  and  18.  In  order  for  a  model  based  on  trapping  by 
bulk  impurities  to  give  the  temperature  dependence  observed  experimentally,  it 
is  necessary  to  introduce  an  energy  barrier  for  electron  capture.  Clearly  a 
model  based  on  Fig.  5  cannot  give  the  required  temperature  dependence.  In 
this  simplest  formulation,  the  capture  rate  of  electrons  is  given,  for 
example,  by 

Re  =  VcmNan  (39) 

which  is  mathematically  equivalent  to  a  statement  that  all  that  is  required 
for  an  electron  and  a  neutral  acceptor  to  recombine  is  that  they  be  in 
sufficiently  close  proximity  to  each  other. 


(38) 


17 


A  model  for  nonradi ati ve  transitions  in  a  crystal  which  can  perhaps 
explain  the  observed  behavior  on  SI  InP  was  first  proposed  by  Huang  and 
Rhys^J  and  expanded  on  by  Kubo  and  Toyozawa^-I  and  by  Ri ckayzen^^^.  In 
this  model  it  is  assumed  that  the  capture  of  an  electron  by  a  deep  level 
impurity  is  accompanied  by  a  local  relaxation  of  the  crystal  lattice.  In  Fig. 
11  the  upper  curve  represents  the  total  potential  energy  of  the  lattice  plus 
free  electron  as  a  function  of  a  generalized  lattice  coordinate  Q.  The  curve 
represents  schematically  a  hyper  surface  in  4-space.  The  lower  curve 
represents  the  total  potential  energy  of  the  lattice  plus  that  of  a  trapped 
electron.  The  minima  of  the  two  curves  are  displaced  due  to  the  lattice 
relaxation  which  occurs  when  the  electron  is  trapped.  There  are  two  processes 
by  which  an  electron  can  transfer  from  the  free  to  the  trapped  state. 
Classically  it  can  undergo  a  transition  if  the  lattice  vibrations  are  large 
enough  to  cause  the  energy  to  cross  point  C  on  the  upper  curve.  Since  the 
system  could  be  in  either  state  and  have  the  same  value  of  energy  and  lattice 
coordinate  there  is  finite  probability  that  the  electron  will  transfer  from 
the  free  to  the  bound  state.  The  activation  energy  for  this  process  will  be 
E0.  At  low  temperature  when  the  transition  probability  via  this  path  becomes 
small,  the  relative  probability  of  direct  tunneling  from  state  Q"  to  A*  and 
the  subsequent  emission  of  one  or  more  phonons  becomes  larger.  At  T  =  0  the 
optical  transition  energy  would  be  EQ.  The  thermal  activation  energy  for 
transition  from  the  bound  to  free  state  is  indicated  by  Ej.  The  applicability 
of  this  model  for  interpreting  measurements  on  GaAs  and  GaSb  has  been 
discussed  by  Henry  and  Langt^l. 

Seemingly  an  immediate  consequence  of  proposing  this  model  for  the 
electron  trap  is  that  one  is  forced  to  abandon  the  concept  that  the  principle 
of  detailed  balance  will  determine  the  relationship  between  the  capture  and 
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emission  rates  at  steady  state  with  no  electric  field  applied.  I  hesitate  to 
use  the  phrase  thermal  equilibrium  to  describe  this  condition  because  it  seems 
as  if  this  type  of  trap  would  not  be  in  thermal  equilibrium  with  the 
conduction  band  electrons.  The  relative  magnitudes  of  the  capture  and 
emission  coefficients  would  determine  the  density  of  occupied  traps  at  steady 
state  and  the  ratio  of  filled  traps  to  the  total  number  of  traps  can  be  used 
to  define  a  quasi -Fermi  level  for  the  electrons  in  the  traps.  That  this 
result  would  be  the  case  can  be  argued  by  imagining  that  the  energy  E0  of  the 
trap  is  increased  as  Ej  is  held  constant.  The  capture  rate  would  then 
decrease  while  the  emission  rate  remained  constant  and  the  density  of  occupied 
traps  would  decrease.  Since  the  Fermi  energy  for  the  electrons  in  conduction 
band  states  would  be  independent  of  the  details  of  the  trap  confi  guration, 
seemingly  the  only  way  out  is  to  describe  the  trap  by  a  different  quasi -Fermi 
level  even  in  the  steady  state  condition.  The  rate  equation  describing  the 
dynamic  behavior  of  this  trap  then  becomes 


dNaj(x) 
— <TE 


e-Eo/kT( 

a 


-Na.(x)K  e^kT 


eaNc 


-etat 

e  Nai(x) 


(40) 


where  the  functional  form  of  the  emission  rate  as  proposed  by  Kubo^S]  has 
been  employed. 

The  calculation  of  the  time  dependence  of  the  channel  current  for  the 
accumulation  MISFET  was  repeated  using  Eq.  35  for  the  dynamics  of  the  electron 
capture  by  the  deep  acceptors  and  Eq.  17  for  the  capture  dynamics  of  the 
shallow  donors.  For  comparison  the  data  of  Lile,  et  al  are  presented  with 
the  results  as  calculated  above  in  Fig.  12.  The  essential  features  of  the 
data  seem  to  be  adequately  represented  by  this  model. 
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The  300°K  current  drift  behavior  of  the  accumulation  MISFET  was 
calculated  using  the  above  model  for  a  series  of  different  deep  acceptor 
concentrati ons  as  shown  in  Fig.  13.  As  can  be  observed  from  this  figure  the 
maximum  permissible  acceptor  concentration  for  current  drift  which  does  not 
exceed  5%  of  the  initial  current  value  is  about  2x10^  cm~8. 

Conclusion 

The  crucial  question  which  can  only  be  answered  by  further  experimental 
work  concerns  the  amount  of  time  required  for  the  deep  acceptor  levels  to  come 
to  equilibrium  after  the  application  of  a  gate  pulse.  If  the  process  occurs 
in  less  than  a  nanosecond  then  the  filling  of  deep  levels  by  electrons  is 
clearly  not  causing  the  drift  in  channel  current  that  has  been  observed.  On 
the  other  hand  if  the  process  takes  several  days  (not  very  plausible)  it  is 
also  not  of  much  technological  interest.  Only  if  the  capture  cross  sections 
are  such  as  to  make  the  effect  observable  over  a  time  interval  from  about  10‘8 
to  10^  sec  will  the  effect  be  important  to  circuit  designers.  In  the  present 
calculations  the  capture  cross  section  has  been  treated  as  a  curve  fitting 
parameter  -  no  independent  measurements  of  capture  cross  sections  for  these 
deep  levels  are  available.  It  is  clearly  important  to  obtain  such  data.  The 
values  of  the  impurity  capture  cross  sections  used  in  these  calculations  are 
considerably  smaller  than  those  quoted  by  Henry  and  Lang^25^  for  GaAs  and 
GaSb.  Until  we  have  such  data  for  InP,  it  will  be  difficult  to  determine  if 
the  values  used  for  the  present  calculations  are  realistic  or  not. 

The  alert  reader  will  have  noted  that  early  in  the  paper  I  made  the  point 
that  the  calculations  of  Okamura  and  Kobayashi  give  the  correct  temperature 
dependence  of  the  current  drift  only  if  they  neglect  the  principle  of  detailed 
balance  in  their  formulation  of  the  rate  equation  to  describe  the  capture  and 
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emission  of  electrons  from  traps  in  the  oxide.  Then  at  a  later  stage  I  have 
invoked  essentially  the  same  mathematical  formalism  to  describe  the  capture  of 
electrons  by  deep  levels  in  the  semiconductor  through  a  multiphonon  emission 
process.  How  do  we  know  where  the  charge  trapping  is  occurring?  Is  it  in  the 
insulator  or  in  the  semiconductor?  Given  the  available  data  the  possibility 
of  drift  due  to  traps  in  the  oxide  cannot  be  eliminated,  however,  it  can  be 
argued  that  at  least  some  of  the  observed  drift  is  due  to  traps  in  the 
semiconductor.  At  very  low  temperatures  (77°K)  we  see  no  drift -on  any  time 
scale.  If  the  semiconductor  traps  were  charging  too  rapidly  to  be  observed  at 
300°K,  we  should  at  least  see  some  evidence  of  their  existence  as  the  sample 
is  cooled.  Instead  we  see  nothing  -  no  drift.  The  calculations  presented  in 
this  paper  have  shown  that  for  deep  acceptor  densities  in  the  range  of  10^ 
cnT'*  (a  lower  limit  for  the  current  materials  technology)  the  effects  of  deep 
level  charging  are  easily  observable.  Therefore  at  least  some  of  the  charging 
effects  that  are  seen  from  77 °K  upward  must  be  due  to  charging  of  states  in 
the  semiconductor. 

In  the  present  calculations  only  the  effects  of  uniform  impurity 
densities  have  been  considered.  If  greater  densities  of  deep  levels  are 
caused  near  the  interface  due  to  processing  induced  damage  then  even  stronger 
effects  would  result  than  those  predicted  from  estimates  of  the  bulk  iron 
density.  Such  results  would  also  be  observed  if  the  processing  were  to  induce 
the  pile-up  of  i ron  near  the  surface. 

Obviously  a  need  exists  for  those  people  working  on  InP  devices  to  begin 
to  carefully  evaluate  the  acceptor  and  donor  densities  in  semi -i nsulati ng 
starting  material  before  trying  to  draw  conclusions  on  the  efficacy  of  various 
surface  preparation  and  insulator  growth  procedures  in  producing  improved  InP 
accumulation-type  transistors.  Such  measurements  are  at  the  moment  time 
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consuming  and  tedious  and  vendors  of  In?  cannot  be  relied  upon  to  provide  such 
information.  However,  there  seems  to  be  no  alternative  if  we  are  to  make 
further  progress  in  developing  this  device. 
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Fi gure  Captions 

Figure  1.  Cross  sectional  view  of  accumulation-type  MISFET  on  InP. 

Figure  2.  Channel  current  as  a  function  of  time  for  an  accumulation -type 

MISFET  on  InP,  after  Li  le,  et  al  .  ^  Channel  length  4  um, 
channel  width  400  um,  oxide  thickness  100  nm,  gate  voltage  4  v, 
drain  voltage  2  v. 

Figure  3a.  Energy  band  diagram  for  the  native  oxide  on  InP  as  proposed  by 
Goodnick,  et  al.E^  with  gate  voltage  equal  to  zero. 

Figure  3b.  Energy  band  diagram  for  the  native  oxides  on  InP  with  a  positive 
gate  voltage. 

Figure  4.  Energy  band  diagram  for  the  tunneling  model  proposed  by  Okainura 
and  Kobayashi 

Figure  5.  Energy  band  diagram  for  an  accumulated  surface  on  semi- 
insulating  InP. 

Figure  6.  Surface  electric  field  vs.  surface  potential  at  t  =  0  and  t  =  » 

for  InP  with  Nfi  =  5x10  cm-3,  Nd  =  lxlO16  cm-3  and  T  =  300 °K. 

The  dashed  lines  represent  the  boundary  condition, 

V  =>|'+eE/C  =  constant,  for  various  values  of  the  gate 

g  s  s  s  ox  a 

voltage,  Vg.  As  time  increases  from  t  =  0  the  solution  for  the 
system  will  move  along  a  line  of  constant  gate  voltage  from  the 
curve  representing  the  initial  state  to  the  curve  representing  the 
final  state. 

Figure  7.  Surface  electric  field  vs.  surface  potential  at  t  =  0  and  t  =  *° 

for  InP  with  Nfl  =  5xl0^17  cm-3,  Nrf  =  lxlO16,  and  T  =  300°K. 

Figure  8.  Plots  of  electric  field  and  surface  potential  vs.  distance  from 
the  semiconductor  surface  at  t  =  0  and  t  =  •  .  Nd  is 
lxlO^3  cm-3,  Na  is  5x10^®  cm-3,  and  T  is  300°K. 
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Figure  9.  Dependence  of  initial  and  first  values  of  the  charge  in  the 

conducting  channel  of  an  accumulation-type  MISFET  on  the  deep 
acceptor  doping.  The  dashed  lines  represent  initial  values  and 
the  solid  lines  represent  final  values.  The  dielectric  thickness 
is  100  nm,  the  relative  dielectric  constant  is  4,  is  10^  cm~^ 
and  T  is  300 °K. 

Figure  10.  Time  dependence  of  the  channel  current  for  an  accumul ation-type 
MISFET  assuming  a  Shockley-Read  type  of  rate  equation,  Vg  is 
4  v,  the  dielectric  thickness  is  100  nm,  the  relative  dielectric 
constant  is  4,  crn(donor  states)  is  1.25x10"^®  cm^  and 
0^  (acceptor  states)  is  1.25x10“^®  cm^. 

Figure  11.  Confi gurati on  coordinate  diagram  for  a  deep  level  trap.  The 

energy  (U)  of  the  total  system  of  the  lattice  plus  one  electron 
is  plotted  as  a  function  of  a  generalized  lattice  coordinate  (Q) 
and  is  plotted  for  the  case  when  the  electron  is  free  (upper 
curve)  and  for  the  case  when  the  electron  is  trapped  (lower 
curve).  EQ  represents  the  energy  barrier  for  electron  capture 
and  Ej  represents  the  energy  for  electron  emission.  The  energy 
for  optical  emission  from  the  minimum  energy.  O',  of  the  system 
when  the  electron  is  trapped  to  the  conduction  band  is  given  by 

e2. 

Figure  12.  Charge  in  the  conducting  channel  of  an  accumulation-type  MISFET 
as  a  function  of  time  assuming  that  electron  capture  by  the  deep 
acceptor  takes  place  via  a  multi  phonon  emission  process.  The 
crosses  represent  data  taken  by  Lile,  et_a]_J-7-^  Nj  is 
l.lxlO*7  cm  ”3,  a^donor)  is  1. 25x10"^  cm^,  cr^  (acceptor) 
is  1.25x10"^  cm^,  and  =  1230  cm^  v”*  s“*. 
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Figure  13.  Charge  in  the  conducting  channel  of  an  accumulation-type  MISFET 
as  a  function  of  time  for  deep  acceptor  concentrations  between 
1015  cm-3  and  1018  cm"3  assuming  that  the  electron  capture  takes 
place  via  a  multi  phonon  emission  process.  Nd  is  1015  cm-3  and 
T  is  300 °K. 


27 


rom  Surface  (meter) 


IEEE  ELECTRON  DEVICE  LETTERS.  VOL.  ED  LA.  NO.  7.  JULY  IMS 


Space-CIianrge-Limited  Currents  and  Trapping  in 

Semi-Insulating  InP 
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Abstract — The  energy  level  and  the  density  of  the  deep  acceptor  Fe3* 
trap  involved  in  space-charge-limited  current  (SCLC)  flowing  in 
semi-insulating  (SI)  InP  was  deduced  from  l-V  measurements  and  a 
model  employing  either  the  low-field  Hall  mobility  or  the  high-field 
saturated  electron  velocity  o,.  The  modified  model  based  on  v,  yields  the 
Fe3*  energy  level,  E,  =  0.6  eV,  and  the  trap  density  N,  =  3  x  10‘Vcm3, 
consistent  with  data  obtained  by  other  means. 

nPHE  INTERACTION  and  parasitic  coupling  between 
-**-  adjacent  field-effect  transistors  on  their  common 
semi-insulating  (SI)  GaAs  substrates  has  led  to  detailed 
investigations  of  charge  carrier  transport  in  these  substrates. 
Measured  current  versus  voltage  ( l-V )  data  have  been 
interpreted  in  terms  of  models  which  include  [l]-[5] 
space-charge-limited  current  (SCLC)  mediated  by  trapping, 
space-charge  domain  formation,  imcact  ionization  in  high 


electric  fields,  and.  more  recently,  self  oscillation  [6],  [7] 
associated  with  field-dependent  trapping.  In  SI  GaAs  the 
main  electron  trap  EL  2  which  is  a  deep  donor  level  and 
compensates  shallow  residual  acceptors  has  a  significant 
effect  on  its  l-V  characteristics.  In  SI  InP  residual  shallow 
donor  levels  are  compensated  by  deliberately  introduced 
Fe3’  deep  acceptors.  Although  qualitative  similarities  are 
expected  between  SI  GaAs  and  SI  InP,  the  differences 
between  the  type  ar.d  characteristics  of  their  deep-level  traps 
were  expected  to  lead  to  observable  differences  in  their 
charge  carrier  transport  properties.  In  view  of  the  applica¬ 
tions  of  SI  InP  for  electronic  aad  electroopric  integrated 
circuit  substrates  investigations  similar  to  those  performed 
on  SI  GaAs  were  undertaken. 

Au-Ge  eutectic  contacts,  350  /an  X  130  /im,  were 
vacuum  deposited  and  annealed  on  the  surface  of 
liquid-encapsulated  Czochralrki-grown  polished 
(100) -oriented  Fe-doped  SI  InP  which  is  slightly  n-type  with 
resistivity  p  =  6  x  107  Q-cm  and  electron  mobility  p  = 
1500  cm2/V's  determined  from  low-field  Hall  measure¬ 
ments.  The  dc  current-voltage  characteristics  were  measured 
at  room  temperature  under  control  of  an  HP  87XM  computer 
using  an  HP  6034A  programmable  dc  power  source  and  a 
Keithley  619  programmable  electrometer/multimeter  to  mea¬ 
sure  the  current.  Fig.  1  shows  the  /-k'data  of  a  representa¬ 
tive  specimen  with  5-pm  separation  between  the  contacts.  It 
is  qualitatively  similar  to  that  of  SI  GaAs  without  the  sharp 
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LOG  V,  volts 

Fig.  I.  Current-voltage  curve  for  Fe-doped  SI  InP,  with  Au-Ge  conucti 
separated  by  S  /un,  showing  characteristic  SCL  form. 


inflection  point  of  the  current  at  the  onset  of  the  trap  filled 
limit  [8J.  The  onset  of  the  trap  filled  limit  is  an  order  of 
magnitude  higher  than  that  obtained  on  SI  GaAs  with 
comparable  contact  spacing. 

Interpretation  of  SCLC  in  SI  GaAs  has  been  based,  as  a 
rule,  on  the  one-dimensional  single  carrier  transport  model 
of  Ln.mpeit  and  Mark  [8]  which  requires  the  assumption  oft. 
particular  trap  level  and  of  its  energy  distribution.  A 
straightforward  method  which  provides  the  trap  energy  with 
no  a  priori  assumption  was  described  by  Pfister  [91  and 
further  developed  by  Manfredotri  et  ai.  [1C).  We  have  used 
this  method  to  analyze  the  data  shown  in  Fig.  1.  Starting  with 
Ohm’s  law  and  the  one-dimensional  Poisson’s  equation 

J  =  epnE 

ddE/dx )  =p  (1) 

where  J  is  current  density,  e  the  electron  charge,  p  the  free 
carrier  mobility,  n  the  concentration  of  free  carriers.  E  the 
electric  field,  t  the  dielectric  permittivity  of  rhe  material,  and 
p  ihe  total  charge  density.  Assuming  a  constant  carrier 
mobility  p.  a  spaiiany  homogeneous  trap  distribution,  and  a 
constant  electrical  length  L  between  contacts,  one  obtains  [9) 


L  \d{V/St)']- 1 

"°~7p  L  d(\/.n  J 


where  na  is  the  free  carrier  density  at  the  anode  ( x  =  L )  and 
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V  the  allied  potential.  The  .quasi-Fermi  level  £>  at  the 
anode  is  then  determined  from 

Nc 

Ef=kT  In  —  (3) 

na 

with  Ef  considered  positive  going  downward  from  the 
conduction  band,  and  Nc  the  conduction-band  density  of 
states,  k  Boltzman's  constant,  and  T  the  absolute  tempera¬ 
ture.  The  total  charge  density  at  the  anode  pa  may  be  found 

from 


_e  d  \d(V/J2)' 

P‘~T2  d(l/J )  L  d(\/J)  . 


If  trapped  charge  dominates  then,  with  n,  =  (1  /e)pa 

1  *1, 

edEF~dEF  U 

where  n,  is  the  density  of  trapped  carriers  at  the  anode. 
(Trapped  charge  dominance  corresponds  to  e~]pa  —  n,  > 
n„).  If  a  discrete  electron  trap  of  energy  E,  and  total 
concentration  N,  is  present  then 


,  1  (EF-E,\ 

,+rn  -tr) 


where  g  is  the  spin  degeneracy  factor.  dn,/dEF  (not  n,\)  will 
have  its  maximum  at 

EFmu  -E,  +  kT  In  g~  E,  (7) 

with  a  value  of 


*/"l  _  N‘ 

4  kT' 


A  simplified  version  of  this  method  [10J  assumes  that  the 
electric  field  at  the  anode  is  essentially  twice  the  “ohmic” 
electric  field  V/L.  In  that  case,  (2)  and  (4)  simplify  to 


respectively.  Using  (3),  (5),  (II),  and  (12)  dn,/dEF  be 
W  expressed  directly  in  terms  of  the  current  density  and  voltage 


dn,  2e 
dEF~ekTL2 


dJ 

J-  V - 

dV 


(Assuming  that  g  =  2  and  that  kT  In  g  is  roughly  0.02  eV  at 
room  temperature  then  it  is  usually  sufficient  to  state  that  E, 
—  EFm„).  Thus,  from  the  position  of  the  maximum  EFm„, 
and  from  the  value  of  the  maximum  itself,  E,  and  N,  may  be 
obtained,  respectively.  A  simple  check  that  Er  is  discrete 
made  by  means  of  (5)  and  (6)  is 


Hence,  if  N,  obtained  from  (8)  agrees  with  that  obtained  from 
(9),  then  the  level  is  likely  to  be  discrete.  An  additional  check 
is  that  the  full  width  at  half  maximum  (FWHM)  of  dn,/dEF 
for  (5)  is 

FWHM  of  dn,/dEF=  3.5  kT  (10) 

which  is  approximately  0.09  cV  at  room  temperature. 

One  drawback  to  the  use  of  >2),  (4),  and  (5)  is  that  multiple 
derivatives  of  the  experimental  l-V  curve  must  be  deter¬ 
mined;  the  errors  and  fluctuations  in  the  data  arc  multiplied 
with  each  derivative. 


This  greatly  simplifies  calculation,  in  that  only  the  first 
derivative  of  the  I-  V  curve  need  be  obtained.  Equations  (3) 
and  (5>— (10)  are  unchanged. 

This  simplified  version  was  applied  to  the  I-  V data,  shown 
in  Fig.  1 ,  with  p  ~  1500  cm’/K-s,  e  =  1.1  x  10~l2F/cm 
(corresponding  to  a  dielectric  constant  of  12.4).  Nc  ~  5.4  x 
lO'Vcm3,  and  L  =  5  pm.  7  was  replaced  by  I/A  where  A  is 
the  contact  area,  350  nm  x  130  pm.  ( dl/dV )  was  calculated 
numerically,  and  some  smoothing  of  the  I-V  data  was 
required  to  attenuate  fluctuations.  Fig.  2  shows  the  resultant 
dn,/dEF  versus  EF  plot.  From  (7).  a  trap  energy  level  of  E, 
—  0.66  eV  is  obtained,  and  from  (8)  and  the  peak  value  of 
dn,/dEF  (6.8  x  I0l6/cm3-eV),  its  density  N,  is  7  x 
lO’Vcm3.  It  is  in  disagreement  with  that  determined  from  n, 
versus  in  the  same  figure,  and  (8),  N,  =  2.4  x  lO'Vcm3. 
Furthermore,  the  FWHM  of  dn,/dEF  is  0.03  eV,  less  than 
the  0.09  eV  expected  fora  discrete  level.  It  is,  however,  in 
good  agreement  with  the  thermal  activation  energy  deter¬ 
mined  by  Mizuno  et  al.  [12]. 

The  electron  velocity  of  InP  reaches  a  peak  value  for  an 
electric  field,  E  =  11  kV/cm,  which  for  a  5-jrm  contact 
separation  is  -5  V.  The  maximum  voltage  of  Fig.  1, 40  V, 
is  well  beyond  this  saturation  limit  and  a  constant  electron 
velocity  vs,  rather  than  a  constant  mobility  may  be  assumed 
applicable  to  SI  InP  with 

J=env,.  (14) 

With  V-  J  —  dJ/dx  —  0,  and  u,  constant,  n  does  not  depend 
on  x,  and 

p=7*V/L2  (15) 

where  the  potential  V  =  —  j  £  Edx.  Equations  (3)  and  (5) 
are  still  applicable,  and  with  (14)  solved  for  n.  there  follows 
from  (15): 

1  dp  dn,  2d  dV 

~  37T~  wrr=TT-71  *  w?  *  <l6> 

e  dkF  dEF  ekTL2  dl 


0690  0 6fi  0662  0648  0634  0620 


El,  eV 

Fig.  2.  Occupied  (rap  concentration  /»,.  and  change  in  occupied  trap 
concentration  with  Fermi  level  dn,/dEf  versus  Fermi  level  Ef,  for  the 
data  from  Fig.  1  assuming  constant  mobility. 
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Fic.  3.  Occupied  trap  concentration  n,t  and  change  in  trap  concentration 
with  Fermi  level  dn,/dEr  versus  Fermi  level  Ef.  for  the  data  from 
Frg.  1  assuming  constant  carrier  velocity. 

This  modified  version  was  applied  to  the  I-V dat3  in  Fig.  1 
and  the  dn,/dE y  versus  £>  plot  obtained  therefrom  is  shown 
in  Fig.  3.  The  saturated  electron  velocity  t/,,  used  was  1.5  X 
"O'  cm/s,  consistent  with  theoretically  calculated  and  experi¬ 
mentally  verified  values  of  n-types  InP  [11]. 

From  Fig.  3  and  the  appropriate  equations  (7)— (10),  E,  = 
0.60  eV  (below  the  conduction  band  edge)  and  N,  =  3  x 
10,!/cm3. 

These  results  may  be  compared  to  the  activation  energy  of 
FeJi  in  SI  InP  of  0.64  eV  determined  by  Fung  et  al.  [12] 
using  photoconductive  measurements  and  with  that  of  Zeisse 
el  al.  [13]  of  0.53  eV  using  galvanomagnetic  measurements. 
The  iron  concentration  in  SI  InP,  determined  from  secondary 
ion  mass  spectroscopic  measurements  [13],  is  3  X 


10,5-10,A/cm3  (within  a  factor  of  two),  in  fair  agreement 
with  our  data.  The  FWIIM  of  dn,/dEy  =  0.07  cV  is  much 
closer  to  the  expected  0.09  eV.  We  tend  to  favor,  therefore, 
the  constant  electron  velocity  data  of  Fig.  3  which  yields  Et 
=  0.6  eV. 

In  conclusion,  a  direct  method  intended  for  the  evaluation 
of  the  energy  and  density  of  traps  involved  in  SCLC  has  been 
applied,  in  simplified  form,  to  Fe-dopcd  SI  InP,  assuming  a 
constant  electron  mobility  and  a  modified  version  for  the 
possibly  more  realistic  assumption  of  a  saturated  (constant) 
electron  velocity.  The  results  favor  a  single,  discrete  trap 
level  located  -  0.6  eV  below  the  band  edge  and  a  trap  density 
of  ~3  x  1015/cm3  associated  with  the  ionized  Fe3*  deep 
level  in  InP. 
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